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The low-lying n(=3,4,5)d ^1)5/2 states alkaline earth ions are of vital importance in a number of 
different physical applications. The hyperfine structure constants of these states are characterized by 
unusually strong electron correlation effects. Relativistic coupled-cluster theory has been employed 
to carry out ab initio calculations of these constants. The role of the all order core-polarization 
effects was found to be decisive in obtaining good agreement of the results of our calculations with 
accurate measurements. The present work is an apt demonstration of the power of the coupled- 
cluster method to cope with strongly interacting configurations. 

PACS numbers: 



The ^5*1/2 — > ^D^/2 transition frequencies in singly 
ionized alkaline earth atoms are among the leading can- 
didates in the search for an optical frequency standard 
[H Qi Q' They are also of interest from the point of 
view of quantum information processing [j| , parity non- 
conservation arising from the nuclear anapole moment 
[3, and astrophysics Q- A knowledge of the magnetic 
dipole hyperfine structure constant [A) for the ^1)5/2 
states of these systems is essential for determining fre- 
quency standards [g]- Hyperfine interactions in these 
states are strongly influenced by electron correlation and 
it is therefore of crucial importance to use a method that 
can treat accurately the strongly interacting configura- 
tions. The primary objective of this paper is to study 
the peculiar trends in the electron correlation in the hy- 
perfine interaction constants for the ^D^/2 states of the 
singly ionized alkaline earth atoms using the relativistic 
coupled-cluster (RCC) theory. This theory is equivalent 
to all-order many-body perturbation theory [3, Q and it 
has the potential to take into account the interplay of rel- 
ativistic and correlation effects. It has been successfully 
applied to understand the role of electron correlations in 
s and p valence systems [l^, llil [l2- However, the ex- 
cited ^i?5/2 states of the alkaline earth ions are strongly 
affected by the core electrons. Systematic investigations 
of the core-polarization (CP) effects in these systems can 
be carried out to obtain important insights into these 
phenomena. 

The A constant is given by ^3 
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where gi — [^] is the nuclear Lande g-factor with /j,/ 
and /ijv being the nuclear dipole moment Bohr magne- 
ton, respectively. The single-particle expression for the 
T operator is given elsewhere [ia |. 

In RCC, the atomic wavefunction j'l/^for a single va- 
lence (u) open-shell system is given by [3, [lj| 



|*.)-e^{l + 54|$.), 
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where |$„) — at|<I>o), with |<I>o) being the Dirac-Fock 
(DP) state for the closed-shell system. In the singles and 
doubles approximation, we have 



T = Ti+T2 = Y, a+aJl + \Y. a^a+abaatZ (0.3) 

a,p ab,pq 

Sv = 5*1^ + S2v = ^ a+GysP + -^ X! "'t^'t^bavsll, (0.4) 
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where a and a^ are the creation and annihilation oper- 
ators. a,b,c.. and p,q,r.. represent core and virtual or- 
bitals, respectively, t^ and t^ are the cluster amplitudes 
corresponding to the single (Ti) and double (T2) exci- 
tation operators in the closed-shell and sfj and s^l are 
the cluster amplitudes for the single {Siv) and double 
(821)) excitation operators involving the valence electron 
(v). Equations to determine these amplitudes with par- 
tial triple excitation effects are given elsewhere (Por ex- 
ample, see Plllj). 

The hyperfine constants are calculated using the gen- 
eral expression for an operator O as 
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TABLE I: Total number of GTOs and active orbitals consid- 
ered in the DF and RCC wavefunction calculations. 



TABLE IIL Dirac-Fock, all order CP and experimental results 
of the magnetic dipole hyperfine constant (A) of the ^1)5/2 
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where we define the effective operators O 
Ny = {^y\Sle^ e^Svl^v). The above expression is cal- 
culated using the method described by [Ifl llil ■ Contri- 
butions from the normahzation factor has been consid- 
ered by 



Norm. = {'^j,\0\'^v){ 
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We use the Dirac-Coulomb atomic Hamiltonian for the 
present calculations. We employ Gaussian type orbitals 
(GTOs) to construct the DF wavefunction as explained 
in 15]. Finite nuclei with Fermi charge distribution are 
considered as given by Parpia and Mohanty [ig. All 
core-electron effects are considered in the DF and RCC 
wavefunction calculations for all the systems. In table 
n we present the total number of GTOs used for the 
DF and active orbitals used for the RCC calculations in 
different systems. In table m we present the gi values 
that are used to calculate A for various systems. 

The terms involving OS2V in the expansion of Eqn. 10.51 
represents the all order CP effect. One could, therefore, 



FIG. 1: Goldstone diagrams representing all order CP effects 
from OS2V Arrows pointing downwards and upwards repre- 
sent core (a) and virtual (p) electrons, respectively. Arrow 
inside the circle represents valence electron. 



consider 05*2^ as an all order dressed CP effect and di- 
agrams associated with them are given in Fig. ^ In 
order to highlight the spectacular role of the CP effects, 
we present the DF and OS2V contributions along with 
the experimental results of A for different states in table 
mil It is obvious from this table that the DF and exper- 
imental results have opposite signs except for Ca+. The 
contributions from OS2V are larger than the DF values for 
all the ions except Be+ which has only one core-orbital. 
Their inclusion significantly improves the agreement of 
A with the experimental values for all the ions. 

A deeper understanding of these results can be ob- 
tained from table Hvl where we have given the all order 
CP contribution to A from different core (a) and virtual 
(p) orbitals. It is clear from this table that the valence 
or the low-lying continuum orbitals and the outer core 
orbitals contribute more than the other orbitals. The 
significant contributions due to the valence orbitals arise 
from the large size of the hyperfine matrix elements and 
the small energy differences involving the valence and the 
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TABLE IV: CP contributions from single particle orbitals. 
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FIG. 2: (color online) Core-polarization contributions to A/gj 
for different systems. 
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FIG. 3; (color online) Arcc/^ values for different systems, 
where the RCC terms correspond to Eqn. 10.51 O (1), OSiv + 
cc (2), OS2.+CC (3), SlOSi. (4), SlOSi. + cc (5), ^^OSa^ 
(6) and normalization (7) terms are taken on X- axis. 



outer core orbitals. Even though the high-lying (s and p) 
continuum orbitals have a large density at the nucleus, 
their contributions are small. Their energies are rather 
large compare to the core orbitals and this makes their 
contributions to A almost insignificant. 

From Eqn. 10.11 it is clear that A depends on both the 
nuclear Lande g- factor (gi) and the matrix element of the 
T operator which depends on the electron wavefunctions. 
It would, therefore, be appropriate to study the variation 
in electron correlation in different systems by examining 
the results of A/gj. 

The CP effects arising from different systems for the 
low-lying ^D^/2 states are presented in Fig. El by plot- 
ting the ratio of the contributions of OS2V with the cor- 
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responding gj. These contributions have negative signs 
and they are very pronounced in heavier systems. 

In tabled we give the results of our RCC calculations 
for the final values of A as well as those of other calcula- 
tions and measurements wherever available. In order to 
appreciate the relative importance of the different phys- 
ical effects, we plot their fractional contributions in Fig. 
Ofor all the ions. The role played by the different physi- 
cal effects is evident from this figure. The overall trends 
exhibited by A for the ^D<^/2 states of the alkaline earth 
ions are indeed very peculiar compare to single d- valence 
neutral atoms [ill and also single s and p valence atomic 
systems [l^, [i3|- This peculiarity can undoubtedly be 
attributed to the overwhelming contributions of the CP 
effects to the hyperfine constants of the states that we 



TABLE V: The magnetic dipole hyperfine structure constants 
(A) in MHz. 



System 


State 


This work 


Others 


Experiment 


SBe+ 


3d ^Dg/a 


-0. 976 






^^Mg+ 


3d "D5/2 


0. 107 


0. 1196 [22j 




«Ca+ 


3d 2D5/2 


-3. 351 


-5. 2 [23] 
-4. 2 [24] 
-3. 552 [25] 
-4. 84 [26] 


-3. 8(6) 


«^Sr+ 


4d ^Dg/a 


2. 156 


1. 1 [27] 

2. 507 [25j 
-2. 77 [2^ 


2. 1743(14) 


"^Ba+ 


5d ^Dg/a 


-11. 717 


9. 39 [26] 


-12. 028(11) 



have investigated. 

Our results are considerable improvements over those 
of previous calculations which were performed by differ- 
ent variants of many-body perturbation theories except 
the latest which are based on relativistic configuration 
interaction (RCI) method |2g. The RCI calculations by 
Itano use a multi-configuration DF extended optimized 



level (MCDF-EOL) single particle basis. He considers 
only a subset of the single and double excitations used 
in our calculations and the core excitations are limited 
only to a few outer core electrons. Also, he has excluded 
the non-linear terms present in our RCC wavefunctions 
given in Eqn. 10.21 For Sr+ and Ba+, the RCI calculations 
give wrong signs for A (see tablelv]!. thereby highlighting 
the extraordinarily strong correlation effects in these two 
systems. 

Our present work on the magnetic dipole hyperfine 
structure constants of the ^1^5/2 states of the alkaline 
earth ions is a testament to the remarkable ability of the 
relativistic coupled-cluster theory to successfully capture 
very strong correlation effects even when other widely 
used methods like the finite order relativistic many-body 
perturbation theory and the relativistic configuration in- 
teraction fail. This feature, if suitably exploited, can 
yield a wealth of very useful information about a wide 
range of the atomic properties. 
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